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Abstract 

A search for the possible precursors of an earthquake and its effect on the data series of diverse 
geophysical parameters has been made in the mountain of Northern Tien Shan. The complex instal- 
lation included the NM64 type neutron supermonitor, detectors of the environmental low-energy 
neutrons, the scintillation gamma-detector, the sensor of the local electric field, a seismograph, and 
a weather-station. The specialized data filtration methodic was elaborated with an enhanced sensi- 
tivity to the transient signals of seismic origin. On the eve of, and after a 5.4 magnitude earthquake 
the fine features in temporal behavior of the intensity of low-energy neutron and gamma-radiation 
background, so as irregularities of the local electric field were observed which seem to be in a close 
correlation with each other. These results may be an evidence of the possibility of experimental 
identification of earthquake's precursors in the time up to 2 — 3 days before the beginning of a 
period of intensive tectonic activity. 

1 Introduction 

The urgency of the short-time forecasting of earthquakes, especially of high-magnitude ones, and 
the search for their possible precursors is caused by the huge social-economic losses resulting from 
disastrous consequences of these natural phenomena. In last time, much attention is paid to inves- 
tigation of the mutual correlations between a number of different geophysical, atmospheric, and 
astrophysical characteristics which may be used for diagnostics of the internal centers of seismic 
tension from the viewpoint of a timely earthquake forecast. 

One of the possible effects being, in particular, under consideration at present time as a probe 
candidate of inner layers of the Earth's crust is the intensity variation of the near-Earth environ- 
mental flux of slow neutrons. It is generally believed that at least a part of this flux, together with 



1 



atmospheric neutrons originating in interaction of energetic cosmic ray particles, constitute the 
neutrons being born through the (<X,n) type reactions of the crust's nuclei with a -particles which 
arise in decays of the natural radioactive elements, primarily radon 222 Rn and thoron 220 Rn. After 
their thermalization in multiple scattering with nuclei, slow neutrons of terrestrial origin can reach 
the surface of the Earth and be detected afterwards [1J. Thus, a perceptible correlation has been 
observed between the local intensity of low-energy neutrons and the Moon phase, which, accord- 
ingly to the most appropriate interpretation, may be connected with a subsidiary effusion of the 
a-radioactive gases in additional stress of the inner crust during the new- and full moon periods 
under the summed gravitational influence of the Moon and Sun [|2l 0. Similar effect on the in- 
ner layers may be caused by accumulated mechanical strain in a seismically active region which 
disturbances, in turn, should be accompanied with corresponding deviation of the environmental 
neutron background from its equilibrium state. Indeed, significant perturbations of the registered 
intensity of neutron counting rate have been repeatedly reported in the period of seismic activity, 
both in local scale fl3l |4j, and globally flU. 

As for the atmospheric neutron component of the cosmic ray origin, in the range of low energies 
its intensity is governed, in particular, by the current state of interplanetary medium. Together 
with direct influence on the intensity of high-energy cosmic ray particles — predecessors of the 
slow environmental neutrons, any abrupt changes in these conditions, like that being observed 
after a solar flare, may trigger the processes of seismic activity inside the stressed regions of the 
Earth's crust which circumstance stipulates another, indirect connection mechanism between the 
earthquakes and neutron variations 

Another geophysical phenomenon which seems to be interesting from the viewpoint of earth- 
quake forecasting is the behavior of atmospheric electricity in preceding period. As it has been 
stated in j6l 13 electromagnetic emissions in a frequency range from kHz to MHz may be pro- 
duced by opening the cracks in creep processes which develop in the mechanically loaded areas 
of the Earth's crust before an earthquake. Possibility to use this phenomenon for prediction of the 
spatial location of a nearest earthquake's epicenter is illustrated by the works |[8ll9l [T0l[TTI . 

The influence of the low-frequency electromagnetic fields on the region of ionosphere — mag- 
netosphere transition may cause anomalous precipitations of the charged particles from the inner 
radiation belt which have been many times registered in a earthquake preceding period from the 
low-orbit satellites |[T2l [T3l [T4l . One could expect that the effects of this kind may cause also 
deviations in the usual behavior of different components of the secondary cosmic rays deep in the 
atmosphere, in particular, in the intensity of the electron- and gamma-radiation background. 

Quite a different approach to the problem of the short-time earthquake forecast is suggested 
in lfT5l . Because of the interrelated and complex nature of seismic events it is advisable to search 
for specific features of correlation pattern between a number of heterogeneous phenomenologi- 
cal parameters (the ground water salinity and the amplitude of geoacoustic emissions in case of 
[15J) which may signal the nearness of an earthquake. The necessity to fix simultaneously a wide 
range of geophysical characteristics for the purpose stipulates a preferable using of the complex 
experimental installations which is just the case of present paper. 

On the date of 1 May, 201 1, a 5.4 magnitude earthquake has occurred in the vicinity of Almaty 
city in Kazakhstan Republic. Accordingly to the data of European-Mediterranean Seismic Center, 
its precise moment is 02:31:29 UT while an aftershock series lasted over the next ten days. The 
effects of this earthquake on a number of geophysical parameters were registered by detectors of 
the Tien Shan Mountain Cosmic Ray Station's complex placed 95.6 km apart from epicenter. In 
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particular, this installation includes a number of neutron detectors, both for the hadronic compo- 
nent of the high energy cosmic rays (the NM64 type neutron supermonitor), and for the low-energy 
neutron background from the nearest environment (detectors of thermal neutrons consisting of a 
set of unshielded proportional neutron counters); the gamma-detector on the basis of a Nal scin- 
tillation crystal which is capable to measure the radiation intensity in the range of 20 — 1200 keV, 
separately in six energy diapasons; the electrometer for registration of the induced charge and 
the ionic conductivity of atmosphere; the high- sensitive detector of the low-frequency (infrasonic) 
variations of atmospheric pressure, and the standard common purpose weather station for registra- 
tion of meteorological parameters (air pressure, temperature, humidity etc). In detail, the system 
of neutron detectors of the Tien Shan mountain station is described in |[T6l and [fTTl ; the system 
of gamma-radiation scintillation detectors — in [18| and l|T9l . The data registration procedure in 
considered measurements consisted of a continuous monitoring of detector signals whereas the 
intensities of neutron pulses, those of gamma- scintillations, and the weather data were stored with 
a 1-min temporal resolution, and the measurements of the electric field — with a resolution of one 
second. 

The data obtained in this seismic event are discussed in the present paper. 

2 Experimental Results and Discussion 

Gamma- radiation. Time history of the intensity of gamma-radiation registered on the stretch of 
±one month around the earthquake is shown in its original form in the upper plot (a) of the fig- 
ure^ The lower energy threshold of detected gamma-quanta in these measurements was 400 keV, 
and the time axis of the plot is centered relatively to the earthquake of 1 May, 2011, its zero 
point coinciding strictly with the earthquake moment. To eliminate accidental interferences of 
electromagnetic nature which sometimes occurred irregularly during the measurements the initial 
experimental data are presented in this plot after being piped through the moving average filter. 

As it is seen in the upper panel of figure [TJ the record of gamma-radiation intensity demon- 
strates two kinds of irregularity. First, there are numerous transitory enhancements with typical 
duration of the order of an hour which are met in multitude on considered plot. This is a com- 
mon effect which usually accompanies precipitation, in particularly rains and hails, and has been 
repeatedly observed before, both in the Tien Shan experiment ll20l and in the measurements of 
other groups (e.g., [|2Tll22l ). Against the end of presented period (second half of May 201 1), these 
enhancements occur more frequently, in connection with the beginning of the summer rain time in 
the mountain of Northern Tien Shan. 

More interesting is another irregularity, a noticeable dip in gamma-radiation intensity seen 
around zero point of the time axis, i.e. in the immediate vicinity of the 1 May, 201 1 earthquake. To 
stress this effect more prominently and to get rid of the short-time fluctuations, the experimentally 
measured intensity distribution was additionally re-calculated using a digital low-pass filter with 
the time constant T = 500 min, and the result of this procedure is presented in the middle (b) plot 
of the figure [T] It is seen both on the plots (a) and (b) that, starting since 29 April, 1.5 days before 
the earthquake, the intensity of gamma-radiation signals has diminished constantly and reached 
its minimum on 30 April, immediately on the eve of the primary seismic shock, while after this 
shock happened an abrupt growth in the counting rate. It should be noted also, that, accordingly to 
accompanying data of the Tien Shan weather station, abundant precipitations have taken place, and 
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the relative humidity remained about 100% during the whole day of 30 April, while on 1-4 May 
there was a period of the fair cloudless weather, so the variation of gamma-ray intensity during the 
days of 29 April-1 May occurred in an opposite direction relatively to its usual behavior in rainy 
days: instead to grow up it drops at a rain and raises in a clear time. 

The bottom (c) plot of the figure [T] demonstrates the continuous record of temporal behavior 
of the intensity of gamma-radiation as it has been regularly registered at Tien Shan station over a 
larger time stretch (of about 0.5 year, since 1 January 2011 and up to the middle of June 2011). 
As before, the initial pulse counts of gamma-detector presented here are smoothed with a low- 
pass filter, and, additionally, the seasonal long-term trend in the data which is generally caused by 
temperature drift of scintillator's photomultiplier tube is subtracted from this distribution. As it is 
seen on the bottom plot, after the filtering procedure being applied the only characteristically sharp 
irregularity in gamma-radiation intensity over the whole presented period remains in the vicinity 
of the earthquake moment (which is put in bold in the figure [jjc)). This uniqueness is an evidence 
against any random coincidence between both rare events. 

Neutron signal. Temporal distribution of the intensity of environmental thermal neutron back- 
ground obtained during the time both preceding and succeeding the earthquake date is shown, with 
a bold curve, in the upper panel (a) of the figure[2j For comparison, the smoothed gamma-radiation 
distribution which has been spoken above is plotted on this panel as well with a thin curve. The 
low-energy neutron detector of the Tien Shan station used in these measurements registered most 
effectively the neutrons in the range of thermal energies, E n ~ 10 -2 — 10 _1 eV. Initial measure- 
ments of the neutron counting rate which have been held with a 1-min periodicity were corrected 
for variation of the atmospheric pressure p accordingly to the formula generally used in the cosmic 
ray variations study [23J: i corr = i ■ exp(0.0012 ■ (p — po)), where i corr and i are the corrected and 
original values of neutron signal intensity, and po = 675 mbar is the mean atmospheric pressure at 
the altitude of Tien Shan station. Similarly to the case of the gamma-radiation signal, the low-pass 
filtering procedure was applied to these data after the pressure correction. 

It is seen in the figure [2ja) that during the 40 — 50 h long period before the 1 May 201 1 earth- 
quake the intensity of low-energy neutrons diminishes which behavior generally repeats the con- 
sidered above behavior of gamma-radiation. (However, any irregularities in this period are totally 
absent in the signal of the high-energy (E n > 1 GeV) neutrons and cosmic ray hadrons which have 
been registered simultaneously by the Tien Shan neutron supermonitor). The mutual likeliness of 
the neutron and gamma-radiation signals immediately before the earthquake moment is illustrated 
by the incut to figure |2fa) which demonstrates a quite high consistency of both distributions with 
correlation coefficient r between them being about 0.94. 

In contrast to the low-frequency time series of the upper panel in figure [2} the distribution plot- 
ted on the figure |2^b) with a thin gray line was obtained from original measurements of the neutron 
counting rates with application of a digital high pass filter calculated over 3000 succeeding data 
points (i.e. having the time constant of 3000 min). This procedure permits to eliminate any slow 
trend from and stress the fast variations in neutron signal. Averaging of this high-frequency dis- 
tribution over 50 neighboring points leads to another envelope curve, shown with a black line in 
figure |2^b) which is imposed on the gray one. As it is seen, the latter curve generally shows an 
obvious diurnal oscillation of neutron intensity; and mostly prominently this oscillation reveals 
itself just in the days which precede and succeed the earthquake moment (which is stressed with a 
bold portion of envelope line in the figure [2jb)). If to take this 3 days long segment (the dates of 
30 April, 1, and 2 May 201 1) and to use it for calculation of an auto-correlation function with the 
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same envelope curve then the maximum correlation coefficient, of the order of 0.65, is achieved, 
again, in the vicinity of the earthquake (see the bold sinusoid-like curve above in the plot of fig- 
ure 2(b)). This analysis shows that some characteristic amplification of the diurnal rhythmic of 
environmental neutron signal has happened in the period of seismic activity. 

Electric field. Another anomaly observed during the days preceding to earthquake was found 
in the behavior of the local electric field which has been measured with a 1 s time periodicity by 
an electrometer facility placed 2.9 km apart from the detector complex of Tien Shan station, and 
91.7 km apart from epicenter. As it is seen in the upper panel of figure[3j the unusual unidirectional 
spike-like outbursts of positive polarity appear some days before the earthquake in the record of 
electric field, simultaneously with the drop of gamma-radiation intensity (shown again with a thin 
curve on the same figure), their frequency growing constantly until the moment of the first shock, 
and disappearing completely at the end of aftershock period. 

The powerful separate enhancements in the strength of electric field, each 3 — 5 s in duration 
(see the plots (b) in figure[3]) all have practically one and the same amplitude, ~six times below that 
which is met usually at lightning discharges in thundercloud. In contrast to usual thunderstorm phe- 
nomena, these anomalous pulses seen in pre-earthquake period have not been accompanied with 
any infrasonic signal from the high- sensitive pressure variation sensor which has been operating 
simultaneously as a part of detector complex. (Remember also that there was a period of a quiet 
clear weather in the day of 1 May, 201 1, without any sign of electric activity in atmosphere). As it 
is seen in the high-resolution incut to the top plot of figure [3} the series of unidirectional positive 
pulses terminates 40 min before the first seismic shock, and reappears only around the moment of 
the most powerful 4.2 magnitude aftershock which has happened on 2:31 UT 2 May, 201 1. Similar 
disturbances of the local electric field observed in a pre-earthquake period have been just reported 
previously (e.g. EUl). 

If to remove all the moments of spike-like enhancement from the temporal record of electric 
field and to subject the remaining data to a low-pass filtering procedure the fine structure of field 
variation may be revealed which is shown in lower plot of the figure [3] Transient irregularities in 
the field variation spectrum are seen on this plot which precede with a nearly one-hour accuracy 
the moments both of the main earthquake shock and a powerful aftershock which has happened 
~6 hours later. 

In comparison of the figures [Tj [2j and [3] one can see that the period of anomalous behavior of 
electric field begins simultaneously with the formerly discussed decrease in the gamma- and ther- 
mal neutron signal, and the most intensive electric activity falls on the moment of the deepest dip in 
the intensity of gamma-radiation and neutron flux, immediately before the main earthquake shock. 
The latter circumstance is another evidence against any random coincidence, and, contrarily, of 
a close correlation between these rarely met anomalies in behavior of the three heterogeneous 
geophysical characteristics. 

3 Conclusion 

Some characteristic influence of a 5.4 magnitude earthquake has been revealed in monitoring 
records of the intensity of environmental gamma-radiation, of the flux of low-energy neutrons, and 
of the local electric field at the complex detector installation of the Tien Shan mountain station. 
Anomalous variations of these quite different geophysical parameters are in a good correlation with 
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each other both during the 2 — 3 days long period preceding the earthquake moment, and immedi- 
ately before and in the time of seismic activity. The results presented here seem to be promising 
from the viewpoint of the forecast problem of powerful earthquakes though further investigations 
in this field still remain to be of a vital necessity. 
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Figure 1: Experimental data: the 1-min resolution time history of the intensity of Ey > 400 keV 
gamma-radiation in its original form (a); same distribution after being smoothed with a 500 min 
low-pass filter (b); the continuous long-term record of gamma-radiation intensity over the whole 
0.5 year long period after smoothing and subtracting its seasonal trend (c). Each time, zero point 
of abscissa axis coincides with the 1 May 2011 earthquake moment. 
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Figure 2: Upper panel: time history of the counting rate of thermal neutrons around the 
1 May 201 1 earthquake moment (bold line) shown in comparison with the gamma-radiation inten- 
sity (thin line); lower panel: diurnal rhythmic and auto-correlation function of neutron signal in 
the vicinity of earthquake (see text). 
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Figure 3: (a): record of the strength of electric field in the time preceding the 1 May 2011 earth- 
quake (thick line) in comparison with the gamma-radiation intensity (thin line); (b): typical struc- 
ture of the short-time field irregularities on the example of two single spikes; (c): fine structure 
of the record of electric field after elimination of the spike-like outbursts and applying a low-pass 
filter. Zero point of the time axis always coincides with the moment of the main earthquake shock. 
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